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We point out and explicitly demonstrate a close connection that exists between featureless Mott 
insulators and fractional quantum Hall liquids. Using magnetic Wannier states as the single-particle 
basis in the lowest Landau level (LLL), we demonstrate that the Hamiltonian of interacting bosons 
in the LLL maps onto a Hamiltonian of a featureless Mott insulator on triangular lattice, formed by 
the magnetic Wannier states. The Hamiltonian is remarkably simple and consists only of short-range 
repulsion and ring-exchange terms. 



I. INTRODUCTION 



The goal of this work is to explicitly demonstrate a 
close connection that exists between two paradigmatic 
strongly correlated systems: a Mott insulator and a frac- 
tional quantum Hall liquid (FQHL). The connection is 
to some degree almost obvious. Fractional quantum Hall 
effect (FQHE) arises when a two-dimensional (2D) liquid 
of interacting charged quantum particles (either fermions 
or bosons), placed in a perpendicular magnetic field, be- 
comes incompressible at certain commensurate filling fac- 
tors, i.e. ratios of the number of particles to the num- 
ber of available degenerate single-particle states in low- 
est Landau level (LLL), which is equal to the number 
of magnetic flux quanta piercing the sample. Mott in- 
sulator is a very similar thing: an incompressible state 
arising at specific filling factors, in this case given by the 
ratio of the number of particles to the number of available 
degenerate localized Wannier orbitals in a given crystal 
lattice. An important difference between a FQHL and a 
Mott insulator is that, while the FQHL is a liquid, i.e. is 
featureless and does not break any symmetries, a Mott 
insulator can be either a liquid or a crystal, i.e. either be 
featureless or break the underlying lattice symmetry. In 
fact, in most cases, at a general fractional filling factor 
(for bosons the filling factor is denned here as the ratio 
of the number of particles to the number of orbitals, for 
electrons it is half that ratio), a Mott insulator will break 
symmetry, as happens, for example, in the parent com- 
pounds of the cuprate superconductors. 1 The connection 
is thus between a FQHL and a featureless Mott insulator. 
This, in our opinion, is the main point that makes this 
connection interesting. Featureless Mott insulators have 
been actively searched for in recent years, both experi- 
mentally and theoretically.- Even though a lot of progress 
has been made, in particular concrete microscopic models 
with featureless Mott insulator ground states have been 
proposed j£ the general ingredients, which are necessary 
in a microscopic model to obtain a featureless Mott in- 
sulator ground state, are not yet known. We believe that 
the FQHL connection may prove to be a useful contribu- 
tion to this field. 

While (at least superficially) rather obvious, the FQHL 
to Mott insulator connection has been largely unex- 
plored. Only very recently it was explicitly pointed out 



in a series of papers by D.-H. Lee et alA and by Bergholtz 
and Karlhede£ It was demonstrated in these works that 
for a quantum Hall system on a torus there exists a 
limit, namely the quasi-one-dimensional (quasi ID) limit, 
reached when one of the dimensions of the torus is made 
comparable to or even smaller than the magnetic length, 
in which the Mott insulator connection becomes simple 
and explicit and the fractional quantum Hall liquid be- 
comes a simple crystal (not a featureless Mott insula- 
tor), with Landau-orbital positions playing the role of 
the "lattice sites" . It was further demonstrated that the 
evolution from the quasi ID to the physical 2D limit is 
(in many cases) smooth, with the 2D fractional quantum 
Hall liquid ground state inheriting the discrete degener- 
acy of the ID crystal, but in the form of a topological 
degeneracy, as the fractional quantum Hall liquid is fea- 
tureless. 

While very elegant and appealing, the picture of 
Refs. [4.5] has some imperfections. The first one is that, 
while the evolution from the quasi ID to the 2D limit may 
be smooth, the 2D thermodynamic limit is still singular 
in the sense that the 2D fractional quantum Hall liquid 
is certainly not a simple crystal with a broken transla- 
tional symmetry that obtains in the quasi ID limit, but 
is a featureless liquid with topological order. Moreover, 
the evolution from the ID to the 2D limits is in fact not 
always smooth: for example it is not smooth in the case 
of the v = 1/2 composite fermion Fermi liquid stated 

It would be more satisfying to have an approach that 
could establish the Mott insulator to FQHL connection 
directly in two dimensions. The main obstacle here is a 
problem with notation. Namely, the standard choices for 
LLL orbital eigenstates, like e.g. Landau-gauge orbitals 
used in Refs. [HQ, are delocalized. This means that the 
energy cost for doubly occupying such orbitals vanishes 
in the 2D thermodynamic limit. Then it becomes hard 
make an analogy to Mott insulator, since the Mott in- 
sulator physics is most easily described in terms of pro- 
hibiting double occupation of some orbitals or nearest- 
neighbor groups of orbitals. This physics is completely 
obscured if one uses delocalized states as single particle 
basis. What is needed to make the FQHL to Mott insula- 
tor connection is a single-particle LLL basis, that would 
consists of functions, localized in all directions in the 2D 
plane, analogous to Wannier functions in insulators. This 
appears to be problematic. It is well-known 6 that con- 
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structing exponentially localized Wannier orbitals in the 
LLL is impossible: exponential localization and a nonva- 
nishing topological invariant, the Chern number, which 
characterizes Landau levels and which is the source of 
the precisely quantized Hall conductance, are incompat- 
ible. For our purposes, however, exponential localization 
is unnecessary: all we need is a basis of normalizable 
orbitals, which have a hnite energy cost of double occu- 
pation. It was explicitly demonstrated by Rashba et alX, 
that it is in fact possible to construct exactly such a basis 
of quasilocalized Wannier-like orbitals, called magnetic 
Wannier functions in Ref. Q : these wavefunctions have 
a gaussian core and a 1 jr 2 tail, 1 /r 2 being the fastest de- 
cay compatible with a nontrivial Chern number. While 
not exponentially localized, these magnetic Wannier or- 
bitals are normalizable and have a finite energy cost of 
double occupation. It will be demonstrated in this pa- 
per that using magnetic Wannier functions as a single- 
particle basis in the LLL, it is possible to map the prob- 
lem of interacting particles in 2D in a strong magnetic 
field onto a problem of interacting particles on triangu- 
lar lattice with one magnetic flux quantum per unit cell, 
described by a short-range time-reversal invariant Hamil- 
tonian. While we believe that the final result is valid, 
with possible minor modifications, for either bosons or 
fcrmions, the arguments, leading to this result, work only 
for bosons. We will thus focus henceforth on the FQHE 
of charged bosons. The above mapping then implies that 
at the filling factors, at which the bosons in the LLL ex- 
hibit FQHE, the ground state of the equivalent model on 
the triangular lattice is a featureless Mott insulator with 
topological order. 

The rest of the paper is organized as follows. In section 
HI1 we review, for reader's convenience, the construction 
of the magnetic Wannier basis and point out why its 
straightforward application to our problem is nontriv- 
ial. In section Mil it is demonstrated that the magnetic 
Wannier basis Hamiltonian has an emergent low-energy 
long-wavelength symmetry, that drastically reduces the 
number of terms in the Hamiltonian and makes it possi- 
ble to construct a simple short-range lattice Hamiltonian, 
faithfully representing bosons in the LLL. In section ITVl 
we explicitly discuss the physical properties of this lat- 
tice Hamiltonian and we conclude with a brief summary 
of the results in section IVl 



II. MAGNETIC WANNIER BASIS 

We will start by reviewing, for reader's convenience, 
the construction of the magnetic Wannier basis, proposed 
in Ref. 7]. One starts from the zero-angular- momentum 
symmetric gauge LLL wavefunction: 



co(r) = 



1 



(1) 



where I is the magnetic length. This wavefunction has 
the form of an atomic-like orbital, centered at the origin. 



To construct a complete set of such atomic-like orbitals 
in the LLL, one can translate the above wavefunction, 
using magnetic translation operators, to sites of any 2D 
Bravais lattice. A priori, the only restriction one can 
place on the form of this lattice is that the unit cell must 
contain exactly one magnetic flux quantum, or in other 
words, its area must be equal to 27r^ 2 . However, it will be 
demonstrated below that in fact the most natural choice 
is the triangular lattice. Translating co(r) to sites of this 
triangular lattice, we obtain: 



c m (r) 
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where ai = ax, a.2 = a(x + v / 3y)/2 are the basis vectors 
of the triangular lattice, m = (mi, m^) with integer mi 
label the lattice sites and Tr = exp[— iR • (p — eA/c)] 
are the magnetic translation operators (we will take the 
charge of the bosons to be — e and assume the symmetric 
gauge A = |B x r). The requirement that the unit cell 
contain exactly one magnetic flux quantum gives |ai x 
&2\ = 27r^ 2 , fixing the lattice constant in our case to be 

a = ^4ir£ 2 /V3. 

The set of functions c m (r) looks similar to a com- 
plete but nonorthogonal set of atomic orbitals in a crys- 
tal. However, this appearance is deceptive, since this set 
of functions in fact possesses a very nontrivial property 
that makes it very different from a simple set of local- 
ized atomic orbitals. This property is embodied in the 
following identity, first established by Perelomov:- 



(-1) 



mi+m 2 , 



,(r)=0. 



(3) 



Eq.© means that the set of functions c m (r) is in fact 
overcomplete by exactly one state. This identity is the 
origin of the nontrivial topological properties of the mag- 
netic Block states, which we will construct below as lin- 
ear combinations of c m (r), namely the nontrivial Chern 
number characterizing the LLL. It will also play a very 
important role in our analysis. 

Given the set of atomic-like wavefunctions c m (r), one 
can follow the standard procedure to construct magnetic 
Wannier functions. One first constructs Bloch functions 
out of linear combinations of c m (r) as: 



*k(r) 



1 



V^V(k) 



^c m (r)e 



ik-r n 



(4) 



Here is the number of degenerate states in the LLL, 
which is equal to the number of magnetic flux quanta 
piercing the sample, and ^(k) is a normalization factor. 
Assuming Bloch functions are normalized to unity over 
the sample area, the normalization factor is given by: 



Kk)=£(-1) 



,/« 2 



(5) 
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The momentum k belongs to the first Brillouin zone (BZ) 
of the triangular lattice and is given by k = fcibj + fc 2 b 2 , 
where bi = (x — y/y/3)/a, h 2 = 2j//av / 3 are the ba- 
sis vectors of the reciprocal lattice. Imposing periodic 
boundary conditions with respect to magnetic transla- 
tions along the basis directions ai i2 , fixes k\^ 2 to be fci.2 = 
2^711,2/ y/N$ with integers m,2 satisfying < 
n i,2 < It follows from Eq.Q that the normal- 

ization factor f(k) vanishes at k = k corresponding to 
(ki, k 2 ) = (tt, 7r). As shown in Ref.[7|, the Bloch function 
at this momentum is still, however, well-defined and can 
be found by carefully taking the limit k — > ko in Eq. (J3]) . 
Magnetic Wannier functions are then obtained from the 
Bloch functions by the inverse Fourier transform: 



4> m (r) 



1 



re 



(6) 



These functions form a complete orthonormal set of 
states by construction. For further in-depth discussion 
of the properties of these wavefunctions see Ref.0. 

Given the complete orthonormal set of magnetic Wan- 
nier functions m (r), we can write down the Hamiltonian 
of interacting bosons, projected to the LLL, using this 
basis. The Hamiltonian has the following general form: 



H = H ( m i m 2l^|m 3 m 4 )6] Ili &] ll2 6 m4 6 m3 , 



(7) 



where creates a boson in a magnetic Wannier state 
<^m( r )> and we will assume the repulsive interaction be- 
tween the bosons V to be a contact interaction: V(r — 
r') = VS(r — r'). The matrix elements in Eq. (J7J can 
be easily evaluated numerically. One finds that all these 
matrix elements are nonzero in the thermodynamic limit 
and are short-range, in the (imprecise) sense of decreas- 
ing in magnitude rapidly with the separation between the 
sites. However, even if one assumes that only the ma- 
trix elements between nearest-neighbor sites may be re- 
tained, one still obtains a very complex Hamiltonian with 
a lot of distinct terms, since the only obvious symmetry 
Eq.([7j) possesses is the symmetry of the triangular lat- 
tice. In its raw form, the magnetic Wannier basis Hamil- 
tonian is then rather useless. It turns out, however, that 
this Hamiltonian does in fact possess a hidden symme- 
try, which is revealed in the low-energy long-wavelength 
limit, in the sense to be defined precisely below. 



III. MAGNETIC WANNIER BASIS 
HAMILTONIAN IN THE LONG- WAVELENGTH 
LIMIT 

To proceed, let us consider our Hamiltonian not in 
the Wannier but in the magnetic Bloch basis, given by 
Eq.Q. Explicitly evaluating the matrix element of the 



contact interaction in the Bloch basis we obtain: 

H = / ( k ' C l' f l , ' f l")' 5 q+q'.q"+Gfcl c+q &k +q ^ k+q »fc k , 
k,q,q',q" 

(8) 

where k, q, q',q" belong to the first BZ, <5 q + q \ q »+G ex- 
presses momentum conservation modulo a reciprocal lat- 
tice vector G and the interaction matrix element is given 
by: 

I(k,q,q',q") 
= V j rfr* k+q (r)vl/* +q ,(r)* k+q „(r)* k (r) 

V/4ir£ 2 



^(k + q)i/(k + q')Kk + q"Mk) 



x e 



r mi,m2,m 3 

-(V« 2 )[rm 1 +^ a +(*m 1 -r m3 ) a -Kr 



i 2 .m 3 / 



x e (i/4^)2-[(r mi +r m2 )xr m3 ] 

X e - i ( k +l)' r '™l e - i ( k + C l')' r 'n2 e i ( k +q")- r m 3 ^ (9) 

We now note the following property of the magnetic 
Bloch functions. Using Eq.Q, we can rewrite Eq.Q for 
the magnetic Bloch function as: 

* k (r) = 1 y c* m (r)e i ^+^ Sxr ^. (10) 



Using the complex conjugate of the Perelomov overcom- 
pleteness identity Eq.©, we then find that the zeros of 
the Bloch function v I'k(r) are located at: 

rmk = r m + - (ai +a 2 ) + £ 2 z x k. (11) 

The zeros of the magnetic Bloch functions thus form a 
triangular lattice, with one magnetic flux quantum per 
unit cell. Different values of the first BZ momentum k 
label different positions of this lattice of zeros relative to 
the lattice formed by the basis magnetic Wannier states. 
Since wavefunctions in the LLL are fully specified, up 
to phase factors, by their zeros, it follows that the Bloch 
functions can be identified with the Abrikosov vortex lat- 
tice states, which form the set of ground states of Eq.® 
at large filling factors. From the viewpoint of the Hamil- 
tonian in the magnetic Bloch basis, the Abrikosov vortex 
lattice states correspond to condensation of the bosons in 
states with a particular momentum k. "Condensation" 
here should be understood in the sense of the Bloch states 
being the solutions of the LLL-projected Gross-Pitaevskii 
equation, which is satisfied by the boson fields 6 k (which 
become c- numbers in the limit of large filling factor): 



dH 

9^-/^ = 0, 



(12) 



where 11 is the chemical potential. The solution of this 
equation, corresponding to the triangular Abrikosov vor- 
tex lattice, is given by: 

/i = 2/(k,0,0,0)|6 k | 2 , (13) 
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which determines the filling factor v = |&k| 2 /A^> in terms 
of the chemical potential. It will be demonstrated be- 
low that 7(k, 0, 0, 0) is independent of k, so that all such 
solutions describe degenerate states at the same filling 
factor, as they should. 

The fact that the magnetic Bloch functions correspond 
to Abrikosov vortex lattice states now leads us to the 
following observation: the functions ^k(r) at different k 
must be related to each other by magnetic translations. 
Indeed, we find the following relation: 



* k (r) = e l7k e^ k r * (r - l 2 z x k) 



(14) 



Here the factor e2 k r is an Aharonov-Bohm phase factor 
from the magnetic translation operator and e l7k is given 
by: 



a '7k _ 



m-? 2 z x k) i v , irit 

*k(0) v/^OMk) m 



x e -(l/40(r m +^2xkr--§k.r m< ( 15 ) 

It is important to note that while 5'k+G( r ) = ^k( r ) as 
it should, e i7k+G ^ e i7k . From Eq.HU) ^ immediately 
follows that the interaction matrix element in Eq. ((8} can 
be written as: 

7(k, q, q', q+q') = 7(0, q, q', q+q')/* (0, q, q')/(k, q, q'), 

(16) 

where all the k-dependence is contained in the function 

/(k,q,q') = e-^^+^+^+o'^^+o+o'-Tk). (17) 

It is clear from the above expressions that all Abrikosov 
lattice states, corresponding to condensation of the 
bosons in states with different k, are degenerate, as they 
should be. 

Let us now see what the Abrikosov vortex lattice states 
correspond to in the Wannier basis. Transforming the 
boson creation operator from the Bloch to the Wannier 
basis: 

bi = -±=Y, b i^ Tm > ( is ) 

V ly 4> k 



one can immediately see that the condensation (in the 
sense defined above) of the bosons into Bloch states cor- 
responds to states with uniform phase winding along the 
basis directions of the triangular lattice in the Wannier 
basis, with the phase gradient given by the momentum k. 
This nature of the Abrikosov vortex lattice states has im- 
portant consequences. 

First consequence, that can be seen immediately, is 
that the imaginary-time action, corresponding to long- 
wavelength boson field phase fluctuations about a given 
Abrikosov state, will lack the usual (V#) 2 term, charac- 
teristic of superfluids, since all states with uniform phase 
gradients have the same energy. Instead, the action will 
have the form (after appropriate rescaling of the time and 
spatial coordinates): 



It then follows that the dispersion of small fluctuations 
around an Abrikosov lattice state is quadratic instead 
of linear (this holds provided the LLL approximation is 



valid) uj 



This result is well-known and has been 



S ~ J drdv [{d T ef + (V 2 0) 2 ] 



(19) 



obtained before by a number of authors^— The above 
derivation of this result, using magnetic Wannier func- 
tions, is probably the simplest and the most physically 
transparent. The fact that the excitation spectrum is 
quadratic, instead of linear, immediately leads one to the 
conclusion^ that Bose condensation or true off-diagonal 
long-range order is absent in this system. This does not 
necessarily mean, however, that the system is not su- 
perfluid: as was shown in Ref.[l2j], the vortices are still 
localized at large filling factors and thus the superfluid 
stiffness is finite. As the filling factor is reduced, however, 
one expects a transition from the Abrikosov vortex lat- 
tice state (a vortex solid) into vortex liquid states, some 
of which will be incompressible quantum Hall liquids. 13 
It is these states that are of primary interest to us. 

The second, and the most important consequence for 
our purposes, is that the lack of the (V0) 2 term in the 
phase action actually follows from an emergent conserva- 
tion law: namely the conservation of the center-of-mass of 
the bosons in any collision process, which becomes exact 
at long wavelengths. To see this, we again return to the 
expression for the interaction matrix element in the Bloch 
basis, Eqs. (|16|) and (|17p . It may seem at first sight that 
all the k-dependent phase factors, that appear in Eq. (fTT|) . 
could be removed by a gauge transformation of the bo- 
son creation-annihilation operators, i.e. &ke* 7k — > &k ? ac- 
companied by the corresponding redefinition of the Bloch 
functions ^k( r ) e_i7k — > v ^k( r )- This is, however, gener- 
ally not possible due to the fact that e J7k does not have 
the same periodicity in the reciprocal space as the Bloch 
functions (the whole function /, of course, does have the 
same periodicity as the Bloch functions). To proceed, we 
will make an approximation: we will assume that we can 
restrict ourselves to configurations of the boson fields, 
corresponding to long- wavelength distortions of the clas- 
sical Abrikosov lattice ground states. This is certainly 
a harmless approximation in the vortex lattice state it- 
self and should remain harmless even for vortex liquids 
as long as the vortex lattice correlation length £ is larger 
than the magnetic length. This is somewhat analogous 
to the semiclassical nonlinear-sigma-model treatment of 
low-dimensional quantum antiferromagnetSfia which can 
successfully describe quantum-disordered states in these 
systems. This approximation implies smallness of the ex- 
citation momenta q, q' compared to the reciprocal lattice 
momenta. We introduce a cutoff scale A for the momenta 
q, q', so that |q|, |q'| < A, and assume that A satisfies the 
inequality l/(<A<l/f. To leading order in the small 
parameter A£ we can then set /(k, q, q') s» 1 and the 
interaction matrix element in Eq. {[8j) becomes indepen- 
dent of k. Transforming the Hamiltonian to the mag- 
netic Wannier basis, we obtain Eq.©, with m labeling 
the magnetic Wannier states and the interaction matrix 
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element given by: 

(mim2|y|m3m4) = gA(mi - m 4 , m 2 - m 4 ) 

+rri2 ,1x13+1x14 7 (20) 

where the function <?a depends on the (very loosely de- 
fined) momentum cutoff A and its explicit form is thus 
rather meaningless (in addition, in a more careful deriva- 
tion this function would be modified by integrating out 
excitations with |q|, |q'| > A). The physically mean- 
ingful information in Eq. ([20)) is contained in the Kroe- 
necker delta-symbol, which expresses the conservation of 
the center-of-mass position of the boson pairs, mentioned 
above, and is a consequence of the approximate indepen- 
dence of the interaction matrix element in (JSj) on k. Such 
a center-of-mass position conservation, but only in one 
spatial direction, is obvious and exact in the Landau- 
gauge orbital basis, where it appears as a direct conse- 
quence of the momentum conservation in the transverse 
direction. 4 In our formulation, first BZ momentum is no 
longer exactly conserved (it is conserved up to a recipro- 
cal lattice vector), since our choice of the single-particle 
basis explicitly breaks translational symmetry (this is a 
price we have to pay for using spatially-localized single- 
particle states). The center-of-mass conservation then 
becomes an emergent conservation law, which becomes 
exact at long wavelengths. 

IV. SHORT-RANGE RING-EXCHANGE 
MODEL ON TRIANGULAR LATTICE 

The center-of-mass conservation law, derived above, 
drastically reduces the number of terms in the Wan- 
nier basis Hamiltonian. The final approximation we will 
make, the justification for which will be provided below, 
is that we can retain only the shortest-range terms in the 
Wannier Hamiltonian. This is a harmless approximation 
provided the characteristic range of the matrix element 
Eq. (|2H)) . which is of order 1/A, is much smaller than the 
vortex lattice correlation length £. This was precisely the 
assumption we made in the argument leading to Eq. (|20p 
and thus Eq. (|20[) and the above approximation are con- 
sistent with each other. 

We then arrive at the following simple short-range lat- 
tice Hamiltonian on the triangular lattice, that we conjec- 
ture faithfully represents interacting bosons in the LLL: 

H = -K^2bl 11 bl n2 b mi b m3 +U^2nl 1 +^2 V mm >n m n m ,. 

P m mm' 

(21) 

Eq. (|2ip is the main result of our paper. The first term 
in (|21[) is the shortest-range ring-exhangc term on the 
triangular lattice (there are three distinct kinds of pla- 
quettes P, as shown in Fig. 1), which is the shortest- 
range and thus the dominant center-of-mass conserving 
pair hopping term. The second term is on-site repulsion 
term (n m = fcJn^m)- The third term represents longer- 
range repulsion. The relevant range of V mm < depends on 




FIG. 1: Three types of smallest-size 4-site plaquettes P in 
Eq. (|21[l on the triangular lattice. Ring-exchange term hops a 
pair of bosons on sites 1 and 2 to sites 3 and 4 and back. 

the boson Landau level filling factor v (i.e. should be 
at least of the order of the mean interparticle distance 
for a given filling factor), and can be restricted to only 
nearest- neighbor repulsive interactions at v = 1/2. 

The sign of K is important and can be fixed by requir- 
ing that Eq. (|2"Tj) reproduce the correct ground state at 
large filling factors, i.e. the Abrikosov vortex lattice. It 
is easy to see that at large filling factors the ground state 
of (|21l) with K > is in fact the Abrikosov vortex lat- 
tice. Indeed, in this classical limit we may replace boson 
operators by c-numbers. The dominant repulsive inter- 
action term in Eq. ([21[) is the on-site repulsion term. In 
the classical limit this will favor equal boson density on 
all the lattice sites. The nature of the ground state will 
then be determined by the ring-exchange term, which is 
the only term in (|2"Tj) . that depends on the phases of the 
bosons and has the form: 

H = -2Kj2cos(0 mi + m2 - 9 ma - m J, (22) 
p 

where 9 m is the phase of the boson field 6 m . It is easy 
to show-i£ that the set of ground states of Eq. (|2"2"|) with 
K > corresponds to all possible uniform phase gradi- 
ents along the basis directions a 4 .2 of the triangular lat- 
tice. As already shown in section IIII1 this corresponds 
precisely to Abrikosov vortex lattice states. Since all 
states with uniform phase gradients are degenerate, one 
also obtains the quadratic dispersion for small phase fluc- 
tuations around any of the ground states. The fact that 
Eq. (f2~Tj) correctly reproduces both the ground state and 
the excitation spectrum of the original boson Hamilto- 
nian Eq.0 at large filling factors reassures us that it 
will in fact faithfully represent Eq.([7]) at all filling fac- 
tors, with properly chosen ring-exchange and interaction 
parameters. 

In particular, let us now consider the case of the filling 
factor v = 1/2. For bosons with contact interaction, the 
exact ground state in this case is the v — 1/2 Laughlin 
liquid^ 

y( Zl ,...,z N ) = ll(z z - Z3 fe-^ 2 ^ 2 . (23) 

i<j 

By our conjecture that Eq. (f2"Tj) faithfully represents 
bosons in the LLL, the ground state of this Hamiltonian 
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at filling factor 1/2 is then a featureless Mott insulator 
with topological order. We expect this to be true when 
K < V, where V is the strength of the nearest-neighbor 
repulsive interactions (U is the dominant interaction en- 
ergy scale and can be taken to be large compared to V) . 
Assuming U K, V, the bosons can be taken to be hard- 
core, i.e. with double occupation of any site of the lattice 
prohibited. Using Holstein-Primakoff transformation 17 
between hard-core bosons and spins of magnitude 1/2, 
we can rewrite Eq. (l2"TT) as the following model of inter- 
acting spins-1/2 on the triangular lattice: 

H = -K Y, S+, S+ 2 S mi S m3 + V £ S z m S* m , . (24) 

P (mm') 

The ground state of Eq. (l24l) at K = has extensive de- 
generacy, corresponding to all possible configurations of 
Sf with at most one unsatisfied bond per every triangu- 
lar plaquette of the lattice. 18 ' 19 When this degeneracy 
is lifted by a small two-spin interaction term of the form 
— J(S^Sj + h.c), the ground state is known to be a 
supersolid, 2 ® i.e. a state which has both long-range order 
in the x, y-components of the spin and a finite- wavevector 
ordering of the ^-components. Our mapping between 
Eq.(f24f and the v = 1/2 fractional quantum Hall liq- 
uid means that the ground state of Eq. (|2"4"|) . in contrast, 
is a spin liquid, i.e. a state with a gapped excitation 
spectrum and topology-dependent ground state degener- 
acy, which is the same as a featureless Mott insulator 
in the bosonic language. When K ^> V, we expect the 
ground state to be a compressible liquid with a quadratic 
excitation spectrum^ most likely a superfluid. 

A somewhat subtle issue, that requires special consid- 
eration, is the issue of the ground state degeneracy of the 
featureless Mott insulator ground state of Eq. (l2~4l and 
the nature of its quasiparticle excitations. The ground 
state degeneracy of the v = 1/2 Laughlin liquid on a 
torus is 2-fold and the quasiparticles are anyons of charge 
±e/2j21 Both the 2-fold degeneracy and the anyonic na- 
ture of the quasiparticles depend crucially on the fact 
that the time-reversal symmetry is broken by the perpen- 
dicular magnetic fiek h 21 ' 22 However, the Hamiltonians 
(12T1) and (l24l) are manifestly time-reversal invariant. The 
information about the time-reversal symmetry breaking 
is contained in the Wannier functions <^) m (r), but not 
in the center-of-mass conserving interaction matrix ele- 
ments (mini2|V|ni3in4), which are all real, as can be 
seen by inspection of Eq. @. The matrix elements, 
which do carry the information about the time-reversal 
breaking, are the ones that do not conserve the center-of- 
mass, as these matrix elements are in general complex. 
The simplest kind of such a matrix element, and also the 
one that has the largest magnitude at short distances, is 
the "correlated hopping" -type matrix element with, for 
example, mi = m3, 1112 ^ 1114. It is clear that such a 
matrix element is, in general, complex. It is also very 
easy to see why such matrix elements are irrelevant at 
long distances (but see below): one simply needs to no- 



tice that 

^(mmi|F|mm 2 ) ~ <5 mi ,m 2 - (25) 

m 

Our main assumption is that this irrelevance continues 
to hold even at low filling factors, such as v = 1/2. This 
should be true as long as the correlation length in a given 
state is significantly larger than the magnetic length. The 
only problem with this is that the topological degeneracy 
on a torus of the incompressible liquid ground state of 
Eq. (l24[) . which is time- reversal invariant, has to be equal 
to 4 (assuming the quasiparticle charge is ±e/2, as in 
the Laughlin state) ^ i.e. double the degeneracy of the 
Laughlin liquid. 

The most natural resolution of this apparent para- 
dox seems to be as follows. The set of 4 degenerate 
ground states of Eq.(f24f must consist of 2 pairs of de- 
generate states, each pair corresponding to the Laughlin 
liquid with the magnetic field directed along z or —z, as 
Eq. (l24[) is invariant under time- reversal. It then follows 
that each such pair of states breaks time-reversal symme- 
try spontaneously. The spin liquid ground state of (j2"4")l 
is then a Kalmeyer-Laughlin-type chiral liquid^ which 
spontaneously breaks parity and time-reversal symme- 
try^ The quasiparticle excitations above such a state 
are charge ±e/2 anyons, as in the Laughlin liquid^ The 
role of the complex center-of-mass nonconserving ma- 
trix elements (mim,2|V| 11131114), which explicitly break 
time-reversal symmetry, such as the correlated hopping 
matrix elements mentioned above, is to act as a small 
"symmetry-breaking field" , that lifts the degeneracy be- 
tween the 2 pairs of states, but is otherwise unimportant. 
One then obtains a 2-fold degenerate ground state on a 
torus with anyonic quasiparticle excitations, exactly as 
in the v = 1/2 Laughlin liquid. This scenario is very 
appealing, especially in light of the fact that there are 
so far only two examples of microscopic models in the 
literature, which have been shown to have a chiral spin 
liquid ground state d 6 ! 27 Both these models, however, are 
significantly more complicated than Eq. (l24l) . Our result 
is of course only a conjecture at this point and needs to 
be verified by an explicit numerical simulation. 

V. CONCLUSIONS 

In conclusion, we have derived an explicit mapping be- 
tween the Hamiltonian of interacting bosons in the LLL 
and a time-reversal invariant Hamiltonian of interacting 
bosons on the triangular lattice with one flux quantum 
per unit cell, Eq. (|2"Tl) . At the filling factors, at which the 
bosons in the LLL condense into incompressible quantum 
Hall liquid states (such as v = 1/2), the ground state of 
this lattice Hamiltonian is a featureless Mott insulator 
with topological order and spontaneously broken time- 
reversal symmetry. The ground state degeneracy of the 
featureless Mott insulator state on a torus is thus pre- 
dicted to be equal to twice the ground state degeneracy 
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of the corresponding Laughlin state, i.e. 4 in the case of 
filling factor 1/2. By the same logic, at odd-denominator 
filling factors, such as v = 1/3, the ground states of ([2Tj) 
should be compressible but non-superfluid liquids ( "Bosc 
metals" corresponding to composite fermion Fermi 
liquid ground states of 2D bosons in magnetic field.— All 
these predictions are testable by either quantum Monte- 
Carlo simulations, since Eq. (|2"Tj) does not have a sign 
problem, or by exact diagonalization of (|2~Tj) . While we 
have demonstrated the FQHL to featureless Mott insu- 
lator connection for the case of interacting bosons, we 
believe that our conclusions also hold, with possible mi- 
nor modifications, in the case of interacting fcrmions as 



well, since the physics of the FQHE and of Mott insu- 
lators does not depend significantly on the statistics of 
the particles. 
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